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We investigate the propagation of the spatial soliton in cylindrical strongly nonlocal media by a novel method of 
image beam of light. The effect of the boundary on the soliton acting as the dynamic force for the soliton steering 
is equivalent to the force between the soliton beam and the image beam. The trajectory of the soliton is ana- 
lytically studied which is in good agreement with the experimental results. © 2009 Optical Society of America 
OCIS codes: 190.0190 



Nonlocal spatial solitons are extensively investigated [1- 
16] since the pioneering work of Snyder and Mitlicliell[l]. 
A topic that has captured a rising interest is the interac- 
tion between solitons and the boundaries [4-8], especially 
in lead glass. The group of Morderchai Segev found the 
boundary caused elliptic solitons and vortex solitons in 
lead glass and discussed the influence of the boundary 
force on the soliton trajectory [5, 6]. In liquid crystal A. 
Alberucci et al demonstrated the power-depended soli- 
ton repulsion at the boundary [7, 8]. The nonlocal nature 
of lead glass lies in the thermal optical nonlinearity (Pois- 
son type) which is intrinsically infinite without bound- 
aries [5]. Therefore the behavior of the nonlocal solitons 
in lead glass can be greatly influenced by the remote 
boundary and one can easily obtain solitons output con- 
trolled by asymmetric boundary forces [6]. In addition, 
the problem of light trajectory is generally settled by 
the method of equivalent particle theory and light ray 
equation [6, 17] since the beam width is quite less than 
the media size. The particle behavior of light suggests 
that the soliton can be taken as a point light on the 
cross section of the propagation media. 

In this paper by analogy with the method of images 
in electrostatics, we introduce a method of image beam 
of light to deal with the problem of the boundary effect 
on the soliton trajectory in cylindrical lead glass. The 
boundary force exerted on the soliton can be equivalent 
to the force of a remote image point light to the soliton. 
The experimental data of soliton steering can be fltted 
by analytical solution in good agreement. 

The system we concern is the light-induced thermal 
self-focusing nonlinearity in lead glass. The heat energy 
of the light is slightly absorbed with an absorption co- 
efficient a and diffused with a thermal conductivity tz. 
The cylindrical boundary of the lead glass is thermally 
contacted by a heat sink at a flxed temperature Tq de- 
scribed by Fig. 1(a). A temperature gradient is yielded 
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Fig. 1. (a) Diagrammatic layout of a soliton launched at 
an offset of xq in a circular cross section of the sample 
with radius R. (b) Sketch map of the normalized refrac- 
tive indexes respectively induced by the source beam and 
the image beam. 

whose distribution is governed by Poisson equation [5]: 
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where T(X, F) is the temperature distribution, 
/(X, Y) = \A{X, y)p is the light intensity with A{X, Y) 
being the paraxial beam and R is the radius of the 
cylinder cross section. The heat transfer equation is 
inherently two-dimensional with circular boundary as 
the soliton is invariant in the direction of propagation. 
The laser-induced temperature change AT = T — To 
gives rise to a proportional increase in refractive index [5] 
An = /3AT with the thermo-optical coefficient and 
in reverse the refractive index change has a strong 
impact on the propagation characteristics of the light 
beam, including the focusing and the steering. The 
propagation of an optical beam in lead glass is moded 
by Equation (1) coupled with the equation for the 
paraxial beam Y)\ 
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where k is the wave vector. We rewrite Eq. (1) and (2) 
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in a dimensionless form: 
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where the normahzed coordinates and functions are set- 
ting as: X = X/R,y = Y/R,z = Z/{kR^),Lp = A/Aq 
with Al = nof^/{af3PR'^) and N = PR'^An/no. It is 
noticeable that N and \(f\'^ in Eq. 3(b) have their coun- 
terparts in electrostatics, potential and charges respec- 
tively. 

We solve the light-induced refractive index distribu- 
tion in the view of Green function method. In circular 
domain, the Green function of Poisson equation can be 
deduced by the method of images [18] which is the sum 
potential of the source charge and the image charge: 



G= ^(In^- 



In- 
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where r{x^y)^r' {x' ^y') and Y"{x"^y") are respectively 
the normalized radius vectors of the field point, the 
source point and the image point and r" = 1/r' de- 
scribed in Fig. 1(b). The image charge simulates the in- 
fluence of all the inductive charges on the boundary, so 
the force applied to the source charge by the boundary 
is equivalent to the force between the image charge and 
the source charge. In our present case, when the source 
light has a beam width small enough compared with the 
boundary size, it can be taken as a point light beam. Bor- 
rowing ideas from the method of images in the electro- 
statics we introduce the method of image beam of light. 
By analogy with the the expression of the charge inter- 
action, we directly write the "force" between the source 
and image beams which provides the dynamic force for 
the steering of the source beam: 
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where (xc,0) is the source beam center, (l/xc,0) is the 
image beam center and y has been set zero since the 
problem is ^-symmetric, p = J \^{x^ — Xc^y')\^dx' dy' is 
the normalized light power which is equivalent to the 
charge density in electrostatics. The 1/r law of interac- 
tion is ever predicted by C. Rostschild et al[3]. 

The interaction of the two beams is mediated by the 
light-induced index which is the solution of Eq. 3(c). 
When we make further quantitative analysis, the solu- 
tion is given by: 
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The integration has been expressed in two terms of A^i 
and N2. Ni is symmetric about the beam center (xc,0) 
since Gi is shift invariant and the beam has a symmetric 



profile (in the general case). It represents the source beam 
induced refractive index in the free space who inversely 
serves as the focusing lens for the source beam to form 
solitons. We can obtain the critical power of the soliton 
Pc = 4:7Tnof<i/{a/3k'^WQ) with wq being the beam width 
by expanding A^i to the second order [19]. The effect of 
N2 on the source beam is equivalent to the effect of the 
refractive index induced by the image beam located at 
the image point (l/xc,0) indicated in Fig. 1(b). 

In the case of wq <C R^ G2 is almost unchanged in the 
profile of the light beam, N2 reads 
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The steering trajectory of a light ray is governed by the 
Eikonal equation [9, 17]: 
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Just as our qualitative anticipation in Eq. (5), the accel- 
eration of the source beam center has a form analogous 
to that of the inter-force between the source and the im- 
age charges. We expand the right hand side of Eq. (8) 
with respect to Xc about Xc = under the condition of 
Xc « 1: 

d'^Xc 
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dz^ 27r ^ 

We only provide the solution under the first-order ap- 
proximation since the solutions under the third and fifth- 
order approximations containing Jacobi cn functions are 
too long to be showed: 



x^^^ = Xcocos(w — ^), 

V TT 
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where Xco is the normalized input offset. The normalized 
output coordinate Xq-' , where / is the normalized length 
of the lead glass, is proportional to Xco when the input 
light power maintains constant. 

Fig. 2(a) pictures the theoretical period of the beam 
steering versus the input offset. The period decreases 
when the input offset approaches to the boundary. This is 
understandable since strong boundary effect accelerates 
the oscillation of the beam. Nevertheless the oscillation 
period is too long compared to our glass length even 
in the strong acceleration case. Fig. 2(b) is the beam 
oscillation at the input offset of 0.8R. 

Our experimental setput is similar to that for the work 
of B. Alfassi et al[6]. The sample is a cylindrical heavily 
lead-doped glass sample with the radius R = 1.5 mm 
and length L = 5 cm (Fig. 1(a)). The absorption co- 
efficient a = 0.07 cm~^, the thermo-optical coefficient 
(3 = 14 X 10~^ the refractive index no = 1.9 and 

the thermal conductivity = 0.7 W/(mK). The 50 fim 
FWHM light input is produced by a double frequency 
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(a) (b) 

Fig. 2. (a) The theoretical oscihation period verses the 
input offset Sohd hne: first order approximation. Sohd 
line: first order approximation; dashed line: third order 
approximation; dashdotted line: fifth order approxima- 
tion, (b) The theoretical period of the beam steering 
versus the input offset under the fifth approximation. 

YAG laser (Verdi 5) with the wavelength 532 nm. The 
soliton critical power is measured to be 500 mW which 
is in good agreement with the analytical value [19] under 
the experimental parameters. The data of the steering 
experiment are obtained under the higher light power of 
700 mW. Fig. 4 gives the comparisons between the ex- 
perimental data and the theoretical fitting curves. The 
first approximation solution is a straight line, it can be 
fitted to the experimental data only when the input offset 
is small enough. The higher order approximation gives 
the better agreement with the experimental result. It is 
noticeable that the propagation distance in Fig. 2(b) is 
far longer than the actual sample size. Yet we can not 
obtain large steering via lengthening the glass or obtain 
the output at the desired period, because the absorption 
will result in the changing of the beam width even the 
collapse of the soliton. 

In conclusion, the soliton propagation in cylindrical 
lead glass is studied. A novel method of image beam of 
light is produced by analogy with the method of images 
in electrostatics. The boundary effect on the soliton is 
treated as the inter-force between the soliton beam and 
the image beam located at the image point. The analyt- 
ical solution of the soliton trajectory is in good agree- 
ment with the experimental results. The method of im- 
age beam of light is useful in the treatment of the bound- 
ary problems, including the influence of the boundary 
force on the soliton trajectory and even the interaction 
between solitons via the boundary effect. 
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